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Abstract 

It is shown that 5D Kaluza-Klein theory stabihzed by an external bulk 
scalar field may solve the discrepant laboratory G measurements. This is 
achieved by an effective coupling between gravitation and the geomagnetic 
field. Experimental considerations are also addressed. 

1 Introduction 

Although the methods and techniques have been greatly improved since the late 
nineteenth century, the precision on the measurement of the gravitational constant, 
G, is still the less accurate in comparison with the other fundamental constants 
of nature [1]. Moreover, given the relative uncertainties of most of the individual 
experiments (reaching about 10^^ for the most precise measurements), they show 
an incompatibility which leads to an overall precision of only about 1 part in 10^. 
Thus the current status of the G terrestrial measurements (see [2]) implies either an 
unknown source of errors (not taken into account in the published uncertainties), 
or some new physics [3]. In the latter spirit, many theories which include extradi- 
mensions have been proposed as candidates for the unification of physics. As such, 
they involve a coupling between gravitation and electromagnetism (GE coupling), 
as well as with other gauge fields present. 
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Here we show that the discrepancy between the present results of the G- measurements 
may be understood as a consequence of the GE couphng. Also, this theory predicts 
a variation of the effective fine structure "constant" a with the gravitational field, 
and thus with the cosmological time (see [4]). 

2 Theoretical background 

An argument initially from Landau and Lifshitz [5] may be applied to the pure 
Kaluza-Klein (KK) action (see [6]): the negative sign of the kinetic term of the 
five dimensional (5D) KK internal scalar field, leads to inescapable instability. 
The question to know which of the two conformally related frames (Einstein-Pauli 
frame or Jordan-Fierz frame) is the best remains debated in the literature, each 
frame having its own advantage. The negative kinetic energy density for the 
field, and thus the instability, occur in both frames. In the following, we perform 
calculations in the Jordan-Fierz frame, where, as we show, the discrepant laboratory 
measurements of G find a natural explanation. Also, it is true that the 5D KK theory 
may yield a zero kinetic term (and thus a zero kinetic energy density, which would 
also be unusual in 4D), but this occurs only when the electromagnetic (EM) field is 
identically zero everywhere. This may be relevant for some cosmological solutions, 
but not for our discussion. On the contrary, we argue here that the <l>-field is tightly 
related to the EM field via the hnk between the compactified space of the fifth 
dimension and the U(l) gauge group. 

Since stabilization may be obtained if an external field is present, we assume 
here a version, KK-^, which includes an external bulk scalar field minimally coupled 
to gravity (like the radion in the brane world scenario) . After dimensional reduction 
{a — 0, 1, 2, 3), this bulk field reduces to a four dimensional scalar field '4> — '4>{x°') 
and, in the Jordan-Fierz frame, the low energy effective action takes the form (up 
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to a total divergence) 

+ J ^[ia^V' a> - t/ - (1) 

where A°' is the potential 4-vector of the EM field, F^g ^ daAg — dp the EM 
field strength tensor, U the self-interaction potential of of the symmetry breaking 
type and J its source term. 

The source term of the ■0-field, J, includes the contributions of the ordinary 
matter, of the EM field and of the internal scalar field For each, the couphng 
is defined by a function (temperature dependent, as for the potential U) fx = 
fx{4> , where the subscript X stands for "matter", "EM" and "$". In addition, 
the necessity to recover the Einstein-Maxwell equations in the weak fields limit, 
implies the following conditions: U{v) = and fEuiv, 1) = f matter {"^ i 1) — 0, 
where v denotes the vacuum expectation value (VEV) of the ■?/'-field. 

The contributions of matter and $ are proportional to the traces of their respec- 
tive energy-momentum tensors. A contribution of the form £0 Jem F^p F'^^ accounts 
for the couphng with the EM field. Besides, though the fit to the data involves 
(f , 1) f ^ AttG/c"^, we may infer that {v , l)v is negligibly small at very 
high temperature (e.g., like in the core of the Sun or at the big bang nucleosyn- 
thesis) and even vanishes beyond the critical temperature (say Tc — 6000 K) of the 
potential U = U{'ijj,T) as one gets v = in that case. 

Following Lichnerowicz [7], let us interpret the quantity 

Geff = § (2) 

of the Einstein-Hilbert term, and the factor £oe// = £o[^^ + '^fEM{"^ , *&) ] of the 
Maxwell term respectively as the effective gravitational " constant" and the effective 
vacuum dielectric permittivity. The effective vacuum magnetic permeability reads 



A*oe// = A*o [^^ + 4:fEM{ip , ]~^, SO that the velocity of hght in vacuum remains 
a true universal constant. Both terms depend on the local (for local physics) or 
global (at cosmological scale) value of the KK scalar field assumed to be positive 
defined. 

The least action principle applied to the action (1) yields the generalized Einstein- 
Maxwell equations and the scalar fields equations 

V„V^^ = -.-|^-| (3) 

and 

$ = - ^ £0 i^"'^ + + + l^'^^V' - ^ ( ) <^>- 

(4) 

3 Solutions in presence of a dipolar magnetic field 

Let us study the spatial variation of $ in the weak fields conditions out of the fields' 
source, but in presence of a static dipolar magnetic field, B — V V{r,(p,9). We 
denote r, (p and 9 respectively the radius from the centre, the azimuth angle and the 
colatitudc. Thus, writing $ = ^(r,ip,d), and taking into account that ^{v) = 
(definition of the VEV) , equation (4) simplifies after linearization as 

At = -^[%!i(..l)« + i^l(W)^ (6) 

Since AV = divB = and {VVf = ^ A{V^) - V AV identically, the solution of 
equation (5) above reads merely 



^=l--[^iv,l)v + ^]V . (6) 



For our purpose, it is sufficient to limit the expansion of the scalar poten- 
tial, V , to the terms of the Lcgendre function of degree one (n = 1) and or- 
der one {m — 1). Hence V — {a^/r"^) [g'l cos 9 + gl sin 9 cos(/? -|- h\ sin 9 sin(p], 



where g^, gl and h\ are the relevant Gauss coefficients, a is the Earth's radius and 



M — ^JiQiY + {g\Y + {h\y denotes its magnetic moment. Setting cos</?i = 
-g\/sJ{g\Y + {h\)\ sin(^i = h\/sJ{g\Y + {h\Y and tan A = ^^/^(^ij^T^, 
the solution of equation (5) then reads (and similarly for ip by making the substitu- 



tion %i ^ - ^) 



1 dfEM ^ ^ ;Uo M 2 



where we have set 

a;(6', = cos^ 6' + cot^ A sin^ 6* cos^ ( 9? + ) — cot A sin26' cos ( <^ + (^1 ). (8) 

Thence, one derives the expression of Geff{r, 9, (p) by inserting the solution (7) above 
in relation (2). 

It is worth noticing that the magnetic potential, V, scales as i?r, where B is 
the magnitude of the geomagnetic field at radius r. Indeed, because of this scahng 
effect, the small spatial variations of the geomagnetic field will influence significantly 
the laboratory measurements of G whereas the large local magnetic fields present in 
the laboratory (e. g., the magnetic suspension used to support the balance beam, 
the magnetic damper, etc..) will not. A rough estimate shows that, even using 
a 30 Tesla superconducting magnet, one still needs to gain at least one order of 
magnitude with the most precise G measuring apparatus avalaible yet. Hence, our 
prediction is consistent with the earlier conclusion of Lloyd [11]. 

4 Comparison with laboratory measurements 

There are presently almost 45 results of G measurements published since 1942 [12]. 
Because of the too numerous uncontrolled systematic biases, the mine measurements 
are excluded from the present study (including them will not change our fit because 
of their lack of precision, typically less than 1%). Also, the more discordant lab- 
oratory measurement (high PTB value [8]) is excluded, since it may suffer from a 



systematic effect. The " official " values are presently G = 6.67259 ± 0.00085 lO'^^ 
(CODATA 86, [9]) and G = 6.670 ±0.010 10"^^ (CODATA 2000, [10]) in MKS unit. 
In the following, all the measurements are weighted equally in the fit. Fitting the 
44 data with these values gives respectively xl = 11-128 and xl = 62.498 {xl = 
per degrees of freedom). If we forget the official values and try a best fit, assuming 
an arbitrary constant value of G, we obtain G — 6.6741 10~^^ kg~^ with 
= 2.255. The fit to the same sample of 44 measurements (figure 1), on account 
of the GE coupling, yields (in MKS units) with = 1.669 : 

^ — = ( 0.149929 ± 0.000017 ) - ( 0.0001509 ± 0.0000252 ) x{L, I). (9) 

10 Ge// 

Prom the above fit, one derives both estimates of the true gravitational constant 

G = ( 6.6696 ± 0.0008 ) 10"^^ kg'^ (10) 

and the coupling parameter 

{v,l)v^{ 5.44 ± 0.66 ) 10-^ fm TeV^K (11) 

The latter quantity, expressed in the canonical form (v , l)v — hcM^'^, yields 
a 5D Planck scale M5 ~ 5.9 TeV of the order of the value that is invoked in the 
literature to solve the hierarchy problem. 
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Location [reference] 


Latitude (°) 


Longitude (°) 




Lower rlutt (^iVlDLj [io, ioj 


/I 1 o 

-4:1.2 


1 ^ A n 


p. ^V/io _i_ n nnnv 
0.0 (42 ± U.UUU 1 

D.D/4D ± U.UUiU 


Wuhan (HUST) [17] 


30.6 


106.88 


6.6699 ± 0.0007 


Los Alamos [18] 


35.88 


-106.38 


6.6740 ± 0.0007 


Gaitnersburg (iNBbj [19, 20\ 


38.9 


-77.02 


6.6726 ± 0.0005 
6.6720 ± 0.0041 


TI>^-.,U^-« /ITT \\ foil 

ooulaer (JiLAj [zlj 


4U 


-iOo.2/ 


D.DoM ± u.uuy4 


Gigerwald lake [22, 26\ 


AC: m T 

46.917 


9.4 


6.669 ± 0.005 (at 112 mj 
0.0 /o ± U.UU^ (at oo mj 
6.6700 ± 0.0054 


Zuricn [z4, zoj 


A'7 A 

4:1 A 


o.oo 


D.D/04 ± U.UUUO ± U.UUiO 

6.6749 ± 0.0014 


Budapest [26] 


47.5 


19.07 


6.670 ± 0.008 


Seattle [14] 


47.63 


- 122.33 


6.674215 ± 0.000092 


Sevres (BIPM) [27, 28J 


48.8 


2.13 


6.67559 ± 0.00027 

Cl dQO _i_ n n 1 1 
O.OOO ± U.Uii 


: f 1 

Fribourg [29] 


A Ci Q 

4o.o 


/.lo 


0.0(04 ± U.004o (Uct. o4j 

6.6735 ± 0.0068 (Nov. 84) 

D.D/4U ± U.UUoo (Dec. o4j 

6.6722 ± 0.0051 (Feb. 85) 


Magny-les-Hameaux [30] 


A n 


2 


0.0 (6 ± U.UUo 


Wuppertal 31] 


Oi.z/ 


/ .io 


O.O/OO ± U.UUii ± U.UUZD 


jjraunscJiweig [r ld) [o. oZ\ 




iU.OO 


0. ( i04U ± U.UUUOD 

6.667 ± 0.005 


Moscow [33, 34] 


55.1 


38.85 


6.6729 ± 0.0005 
6.6745 ± 0.0008 


Dye 3, Greenland [35] 


65.19 


-43.82 


6.6726 ± 0.0027 


Lake Brasimone [36] 


43.75 


11.58 


6.688 ± 0.011 



Table 1 : Results of the most precise laboratory measurements of G published 
during the last sixty years and location of the laboratories. 
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Figure 1: Laboratory measurements with relative uncertainty ^^J^ < 10 ^ and 
measuring time At < 200 s (sample SI, 17 points [16], [17] - [20]"' [23], [14], [29] 
- [31], [32] - [35]). The line indicates the best fit Giab versus the mixed variable x 
{xl = 1.327). Assuming a constant G would yield a bad fit to the data {xl = 3.607), 
mostly because of the HUST value. 
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Figure 2: Giab versus x (whole sample plus the PTB 95 value, 45 points [15] - [37]). 
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Sample 


HO 


HI 


SI 

17 points 
(Fig.l) 


xl = 3.607 (best fit) 
Xy = Zl-oZo [mean oi UUUAiA obj 
xl = 141.46 (mean of CODATA 2000) 


X^ = 1.327 


Whole 
[16] - [37] 
44 points 

(Fig.2) 


xl = 2.255 (best fit) 
xl = 11.128 (mean of CODATA 86) 
xl = 62.498 (mean of CODATA 2000) 


xl = 1.669 



Table 2 : Reduced x^ ^^r the two different hypothesis HO (Hypothesis of a con- 
stant G) and HI (Hypothesis of an effective G), and different samples SI and whole 
(except the high PTE value [8], see text). 

Considering the whole sample, we check the relevance of our result under HI 
compared to the best value of G under HO, by applying the F test (Fisher law). 
This yields F^^ = = 16.09, which indicates that, independently of the number of 
parameters (two instead of one), our fit is better with a significance level greater than 
99.9% [13]. Let us emphasize that the most precise value of G today [14] contributes 
to xl to less than 0.0006 in the above fit. Likewise, the last pubhshed value of G 
[38] would contribute to less than 0.02. This suggests that the agreement may be 
better than purely indicated by the values. Moreover, if one substitutes G = 
(6.6731 ± 0.0002) 10"^^ kg'^ s'^ [39] announced by the MSL team in June 2002 
(CPEM, Ottawa, Canada) for the prior G = (6.6742 ± 0.0007) 10"^^ kg-'^ s'^ 
published in 1999 [16], one would obtain with the sample SI : 

m : xl = 5.168 

HI : = 1.162 

and with the whole sample of 44 measurements : 
HO : = 2.836 
m:xl = 1.488 
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5 Discussion and Conclusion 



It is worth noticing that the scalar fields under considerations identify neither to the 
dilaton nor to the inflaton of higher dimensional theories, without further assump- 
tions. In particular, the computation of both scalar fields, as given by equations (3) 
and (4), involves only (sec the right hand sides) quantities related to fields sources, 
and not to the test bodies. Hence, the effective Ge// given by relation (2) does not 
depend on the composition of the test bodies. Besides, let us emphasize that the 
equation of motion of a neutral point-like particle in the genuine KK theory reduces 
to the 4D geodesic equation after dimensional reduction [40], although the KK scalar 
field is coupled to the Maxwell invariant F^pF"'^. In a forthcoming paper [41], we 
address the effect of the varying effective couphng constants on the masses of com- 
posit particles. On account of the Higgs mechanism of quarks and leptons masses 
generation, by promoting the Yukawa coupling constants to effective parameters (on 
an equal footing with G or a) that depend on both scalar fields ip and we prove 
that the KK'^ model is actually consistent with the current experimental bounds 
on the violation of the equivalence principle. Hence, we conclude that present lab- 
oratory experiments may not measure a true constant of gravitation. Instead, in 
addition to all other possible biases (e. g., anelasticity in the wire of torsion pendu- 
lum as pointed out by Kuroda [42] , and which since has been generically taken into 
account in the experiments), they may be pointing out an effective one depending 
on the geomagnetic field at the laboratory position. 
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